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Abstract 

We investigate flow in semi-peripheral nuclear collisions at AGS and SPS energies 
within macroscopic as well as microscopic transport models. The hot and dense zone 
assumes the shape of an ellipsoid which is tilted by an angle @ with respect to the beam 
axis. If matter is close to the softest point of the equation of state, this ellipsoid expands 
predominantly orthogonal to the direction given by 0. This antiflow component is 
responsible for the previously predicted reduction of the directed transverse momentum 
around the softest point of the equation of state. 

Transverse collective flow in relativistic nuclear collisions reveals the properties of the 
nuclear matter equation of state far from the ground state ^ ||, |^, ^, ||, ||, ^, H, ^ . In principle, 
one can distinguish three different types of transverse collective flow: radial, directed, and 



elliptic p. Recent data on directed and elliptic flow |T0|, |Tl], [T2|, |T3[ has revived theoretical 
interest in this subject @, 0, §]. 

Directed flow occurs only in semi-peripheral nuclear collisions, and therefore must be 
studied in fully 3+1-dimensional geometries. The beam axis is in general taken to be the 
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2;-direction, and the reaction plane to be the z — x-plane. At BEVALAC energies, the two 
nuclei "bounce off" each other, giving rise to a positive average momentum {px{y)/N) per 
nucleon in the forward direction [Q. In momentum space, the flow of matter can be described 
in terms of an ellipsoid, defined by the principal axis' of the tensor of inertia [0, |I|, which 
is tilted in the reaction plane by an angle Ofiow with respect to the beam axis. However, the 
actual shape of the distribution of matter in momentum space needs not be ellipsoidal, see 
below. 

In this paper we show that the situation is fundamentally different if the equation of state 
of nuclear matter is softened, either by a phase transition to the quark-gluon plasma or by 
the creation of resonances and string-like excitations. To this end, we employ one- |^ and 
three-fluid dynamics ||^, as well as the microscopic model UrQMD We demonstrate 



that, around AGS energies, the event shape resembles an ellipsoid in coordinate space, tilted 
by an angle G with respect to the beam axis. This ellipsoid expands predominantly ortho- 
gonal to the direction given by B; we therefore term this flow component antiflow. Around 
midrapidity, the antiflow largely cancels the directed flow from the "bounce-off" of the two 
nuclei [0]. We emphasize that here "antiflow" does not mean the flow of antiparticles [p!7[| . 



which is an absorption phenomenon, nor the low energy [i.e. E^^^ ~ 100 MeV/N) antiflow 
due to attractive potentials ||18| . 

This antiflow component has impact on studies of transverse elliptic flow within sim- 
plifled geometrical overlap models ||, |T9[. These studies assume that the longitudinal 
flow vanishes at ^ = in the whole transverse plane. The non-trivial ellipsoidal event shape, 
however, couples longitudinal to transverse flow, and the longitudinal flow no longer vanishes 
everywhere in the transverse plane at z = 0. The amount of longitudinal flow is sensitive to 
the equation of state, as well as the impact parameter and the bombarding energy, and can 
only be determined in fully 3-|-l-dimensional calculations. 

In order to measure the EoS, i.e., in fluid-dynamical terms the pressure p{e, p) as a 
function of energy density e and baryon density p in the local rest frame of a fluid element, 
one studies the transverse momentum in the reaction plane, px- This quantity is proportional 
to the pressure created in the hot and dense collision zone 



J pA^ dt 



The pressure p is exerted over a transverse area A^^. For increasing bombarding energy, the 
flow, ~ Px, flrst increases, as the compression and thus the pressure grow. However, at large 
E^^^ the time span of the collision decreases, diminishing the flow again. The flow is thus 
maximized at some intermediate bombarding energy. 

A phase transition softens the EoS |@]. The pressure increases slower with e and p than 
in the case without phase transition, reducing the velocity of sound. This delays the fluid- 
dynamical expansion considerably, giving the spectators time to pass the hot and dense 
zone, before they are deflected. One-fluid calculations therefore show a local minimum 
(at ~ 8 AGeV) of the excitation function of the directed flow per nucleon, deflned as 
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This is the weighted mean transverse in-plane momentum {px/N{y)) per nucleon, introduced 
in [jl4|. The weight is the net-baryon rapidity distribution, dN/dy. In a fluid-dynamical 
context, the mean transverse momentum {px/N{y)) is defined as 



/Px. A /j,d^xi?(x)mjvM^(x) 

(iv^^V = J^d^^) ' 

Ux = jVx denotes the x-component of the local 4- velocity field, and rriN- is the nucleon 
rest mass. R is the zero-component of the net-baryon 4-current, R = jp. Here, thermal 
smearing is neglected, and it is assumed that the x-component of the nucleon momentum 
can be approximated by m^v Ux- The volume integration is performed over all fiuid elements 
(projectile and target) around a given rapidity y.^ 
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Figure 1: Time-evolution (in the CM-frame) of directed fiow, p'^"/N, for a Au + Au reaction 
at 8 AGeV, 6 = 3 fm, with and without phase transition to the QGP, calculated in one-fiuid 
dynamics. The numbers in parentheses denote the mean net-baryon density in units of the 
ground state density po — 0.16 fm~'^ at the end of the time evolution. 




The EoS used in our one- and three-fiuid calculations includes a first order phase tran- 
sition to a quark-gluon plasma (QGP). The hadronic phase consists of nucleons interacting 
via exchange of a and u mesons and of non- interacting, massive pions. The QGP phase 
is described in the framework of the MIT-Bag model ||2l| as a non-interacting gas of massless 
u and d quarks and gluons, with a bag parameter B^^^ = 235 MeV, resulting in a critical 

^ In the three-fluid model, since the third fluid is by construction baryon-free, the integration covers only 
projectile and target fluids. 
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temperature Tc — 170 MeV. There is a first order phase transition between these phases, 
constructed via Gibbs' conditions of phase coexistence. 

In Fig. |l], we compute the time evolution of the directed flow, p'^" /N, in one-fluid dynam- 
ics, for a Au+Au colhsion at impact parameter 6 = 3 fm and collision energy E^'-^ = 8 AGeV. 
One observes that, due to the softening of the EoS in a phase transition to the QGP, less 
directed flow is produced in the early compression stage than in a purely hadronic scenario. 
In contrast to the hadronic case, where the directed flow remains constant after reaching its 
maximum, in the case of a phase transition, the directed flow decreases again. By the time 
the mean density drops below nuclear ground-state density, pt" /N is reduced to ^ MeV. 
If one follows the fluid evolution even further (to unphysically small values of the density), 
p'^"/N becomes negative. 
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Figure 2: Net-baryon density R (for the same reaction as in Fig. I) at t = 12 fm/c in the 
reaction plane with velocity arrows for midrapidity {\y\ < 0.5) fluid elements: Antiflow - 
thin arrows, Normal flow - bold arrows. 



This observation is explained by an antiflow component which develops when the ex- 
pansion sets in. This phenomenon is shown in Fig. |^, which is a contour plot of the baryon 
density R, with arrows indicating the fluid velocity. Normal flow (bold arrows) is positive 
in the forward hemisphere, and negative in the backward hemisphere, respectively. On the 
other hand, antiflow (thin arrows) is positive in the backward hemisphere, and negative in 
the forward direction. We show velocity arrows for fluid elements within ±0.5 units around 
midrapidity, since this phenomenon develops at midrapidity, as discussed in detail below. 
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Similar results have been reported in P2[ within the microscopic quark-gluon string model 
2|. 



In Fig. one observes that the hot and dense zone assumes the shape of an ellipsoid 
tilted with respect to the beam axis by an angle 0. The ellipsoid expands preferentially 
in the direction where its surface area is largest, cf. ([^), i.e., orthogonal to the direction 
of the normal flow. This causes the antifiow. Moreover, expansion into the direction of 
normal flow is blocked by the spectators. (Similar arguments led to the prediction of in- 
plane eUiptic flow at high bombarding energies ^ ^.) Note that, at z = 0, antifiow has 
a negative longitudinal component for x > 0, and a positive component for x < 0. This is 
the aforementioned coupling of longitudinal and transverse flow in the central plane. 

The evolution of the distribution of nucleons in momentum space is depicted in Fig. ^ 
In this one-fluid calculation, the participants are shifted instantaneously to midrapidity. In 
the early stage, t ~ 2.4 fm/c, they can be found around p ~ 0. (For clarity, the Fermi- 
momentum and the thermal momenta of the nucleons in the local rest-frame of the fluid are 
not included.) At t ~ 4.8 fm/c the normal flow builds up around central rapidities, leading 
to an 'ellipsoidal' distribution. The principal axis is tilted with respect to the rapidity 
axis by Oaow — 7r/4. However, at even later times, an additional orthogonal component, 
the anti-flow, builds up. This is due to the expansion of the above-mentioned ellipsoid in 
coordinate space, which proceeds in the direction of maximal surface. The final distribution 
in momentum space can even be dominated by the anti-flow component and therefore does 
not exhibit an ellipsoidal shape. 

To illustrate the three-dimensional structure of the expanding matter in coordinate space, 
we also show the baryon density distribution in the transverse plane at various values of z. 
The system in the central plane, at z = (upper left panel) is symmetric around the reaction 
plane, confirming the assumption made in simple geometrical overlap models used to study 
elliptic flow (however, as discussed above, the longitudinal flow does not vanish at ^ = 0). 
Furthermore, for z > (forward, or projectile, hemisphere) the system is displaced towards 
positive X. For z = 1.275 and 2.475 fm, antifiow is clearly visible as flow of matter towards 
negative x. In the most forward plane, z = 6.075 fm, only spectators remain, which are "cut 
off" from the central region and flow mainly in the positive x-direction. 

We have also studied directed flow in the three-fluid model, with a dynamical local uni- 
fication procedure. The three-fluid model treats the nucleons of the projectile and target 
nuclei as two different fluids, since they populate different rapidity regions in the beginning 
of the reaction. The same holds for the newly produced particles around midrapidity, which 
are therefore collected in the third fluid. Thus, the three-fluid model accounts for the non- 
equilibrium situation during the compression stage of heavy-ion collisions. The coupling 
between the projectile and target fluids leads to a gradual deceleration and is parametrized 



by free binary A^A^-collisions p4 |. 



The unification of fluids i and j consists of adding their energy-momentum tensors and 
net-baryon currents in the respective cells, 

rrix) + T^ix) = TZUx) , Kix) + Nf{x) = iVLfied(^) (4) 
and common values for e, p, p and are obtained from T[^Z^^ = (e + p)u^u^ — pg^'^, 
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Figure 3: Evolution of the net-baryon number in momentum space within one-fluid dynamics 
(for the same reaction as in Fig. p. 
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^unified = pu^, and the given EoS p = p{e,p). The local criterion for unification is 

P^±Pl>0.9 . (5) 
P 

Here, pij denotes the pressure in T-'J , and p the pressure in T^j^jgg^. Eq. (^) has already been 
used in as a measure for the equilibration process. 
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Figure 5: The excitation function of directed fiow p'^^ /N for Au + Au collisions at impact 
parameter 6 = 3 fm. Dotted lines are results from a one-fiuid calculation; triangles are 
for a purely hadronic EoS, circles are for an EoS with phase transition. Solid lines are 
calculated with the three-fluid model, with (large circles) or without (small circles) dynamical 
uniflcation. All three-fluid calculations are performed with an EoS with phase transition. 



Fig. ^ shows the excitation function of directed flow p'^"/N calculated in the three-fluid 
model in comparison to that obtained in a one- fluid calculation [Q]. Due to non-equilibrium 
effects in the early stage of the reaction, which delay the build-up of transverse pressure 
1^, ^, the flow in the three-fluid model is reduced as compared to the one-fluid calculation 
in the AGS energy range. Furthermore, the minimum in the excitation function of the 
directed flow shifts to higher bombarding energies. The case without dynamical uniflcation 
yields the least amount of stopping and energy deposition, while the one-fluid calculation 
has instantaneous full stopping and maximum energy deposition. The three-fluid model 
with dynamical uniflcation lies between these two limits; it accounts for the limited stopping 
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power of nuclear matter in the early stages of the collision and mutual equilibration of the 
different fluids in the later stages. Consequently, the shift of the minimum is large without, 
and rather moderate with unification. 

The three-fluid model predicts a local minimum in the excitation function of directed 
flow at bombarding energies between 10 and 20 AGeV, depending on the fluid uniflcation 
criterion (j^). While measurements of flow at AGS energies [0 have found a decrease of 
directed flow with increasing bombarding energy, a minimum has so far not been observed. 
In the three-fluid model with unification, the directed flow exhibits a local maximum at 
-^Lab ~ 40 AGeV. If recent CERN-SPS experiments find larger values for the directed 



flow than at the maximum AGS energy, the existence of a minimum in the excitation function 
of the directed flow due to the intermediate softening of the EoS would be unambigously 
proven. 

We emphasize that the excitation function depicted in Fig. ^ has been calculated for fixed 
impact parameter b, which is not directly measurable in an experiment. Usually the amount 
of transverse energy or the number of participating nucleons are employed as measures for 
b, assuming that the interaction volume is given by the geometrical overlap of two spheres 
displaced in x-direction by the amount b. However, the above discussion suggests that such 
a geometry is oversimplified. The two nuclei are partly defiected and stopped, and thus do 
not penetrate as deeply as compared to the simple geometrical overlap case. Furthermore, it 
is also not obvious that the same Et / E^'^^-hin at different bombarding energies corresponds 
to the same impact parameter, since the system geometry may change considerably due to 
energy-dependent phenomena like the stopping power, phase transitions, etc. This should be 
kept in mind when considering the different values of b where the directed flow is strongest: 
6m ~ 4 fm at AGS and 6^ ~ 8 fm at SPS |Tl|. A detailed study of the impact para- 



meter dependence of directed flow, transverse energy production and number of participating 
nucleons within the three-fluid model is in preparation. 

Let us return to the discussion of the antiflow, which develops also in the three-fluid 
model at energies around the minimum in the excitation function of p'^}^ /N. It leads to a 
plateau in {px/N{y)) around midrapidity. In Fig. ^, this quantity is shown as a function 
of rapidity y at different times. Observe that in the late stage of the reaction, close to 
freeze-out, the flow around y = is even negative. 

We locally decompose the flow into a normal component and an antiflow component 



normal flow : |/(x)pa;(x) > , (6) 
antiflow : ?/(x)p^.(x) < . (7) 



Consequently, we deflne 



/pflow/antiflow \ / d^X i?(x) M, (x) [±y (x) (x)] 

The individual components {vl^^/N^y)) and (p^"*'^°^/A^(?/)) are also shown in Fig. ^. The 
antiflow component develops from midrapidity after ~ 6 fm/c. According to the deflni- 
tion (Rl), both the normal flow and the antiflow are discontinuous at ?/ = 0. The sum of 
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Figure 6: Mean in-plane momentum per nucleon {Px/N{y)) as function of rapidity y at 
various times (dots). The dotted and dashed hues show the decomposition into flow and 
antiflow, respectively. 



10 



both components yields the total flow, eq. (0), which is continuous (and equal to zero) at 
midrapidity. 

This phenomenon is not only limited to fluid-dynamical models. The plateau around 
midrapidity in {px/N{y)) is also visible in the microscopic UrQMD model [0. Fig. |^ 
shows the respective {Px/N{y)) for Au + Au collisions at various bombarding energies. The 
flattening around midrapidity is more pronounced at larger energy. A similar behavior has 
been found in other microscopic models p2| , Unlike the fluid-dynamical calculations, 
the microscopic transport model does not really show a negative slope of {Px/N{y)) around 
midrapidity. Unfortunately, measurements of {px/N{y)) at beam energies E = 6 — ll AGeV 
cover only the range y > 0.3 yp^cMS |^ ill- The quantitative value of antiflow at midrapidity 
thus remains undetermined. 
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Figure 7: In-plane transverse momentum distributions at various energies as calculated in 
the UrQMD model. Calculations were performed using density-dependent Skyrme potentials 
corresponding to a hard EoS with an incompressibility K ^ 380 MeV. 



Besides the transverse directed momentum pf^ /N, many other measures for the in- 
plane flow exist, for instance, the slope of {px/N{y)) or {px/N{y/y\,ca.m)) at midrapidity, or 
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Figure 8: Rapidity distributions of normal flow (fuU line) and antiflow (dashed line) at 
various times. 
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their maximum values as function of ?/, or their values at ?/beam/2. All these measures will 
yield qualitatively and quantitatively different results, because flow and antifiow develop in 
different rapidity regions, as seen in Fig. |^. In the early compression stage of the reaction 
the spectators near projectile and target rapidity are deflected by the pressure in the central 
hot and dense zone, producing normal flow (Fig. full line). When the expansion of the hot 
and dense zone sets in, both normal flow and antifiow (Fig. ^, dashed line) develop around 
midrapidity. Nevertheless, the antifiow finally occupies a broader region around midrapidity 
than the normal flow, while the normal flow dominates the region near projectile and target 
rapidities. 

In summary, we investigated transverse directed flow in macroscopic as well as micro- 
scopic transport models. For the three-fluid model, we find a minimum in the excitation 
function of the directed flow at bombarding energies between 10 and 20 AGeV, which is 
somewhat above the value found earlier in one-fluid calculations. The minimum is caused 
by a softening of the EoS due to a phase transition to the QGP. An antifiow component 
was identified as source for the reduction of the directed flow and discussed in detail. We 
also found that the directed flow of nucleons increases again at higher bombarding energy, 
leading to a maximum in the excitation function at energies around 40 AGeV. If the directed 
flow at this energy, currently investigated by CERN-SPS experiments, proves to be larger 
than at maximum AGS energies, the existence of a minimum in the excitation function of 
the directed flow would be unambigously proven. 
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